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O ■ Abstract 
(N 



^ ■ Channel side information (CSI) information varies very fast in the milhmeter wave radio frequencies so 

its coherent time is much smaller than that of the CSI in low-frequencies (about 6 GHz or below) under the 
same mobility conditions. Further, when the feedback channels, regardless of being wireless or wireline, have a 
backhaul /transmission delay which is practically so much higher than the coherent time of the CSI in the millimeter 
wave radio frequencies. As such in practical scenarios, the CSI will be inherently changed in the fading environment 
before successfully reaches to the central entity (CE) through feedback channels, which specifically means that the 



^SJ ■ CSI reached at the CE certainly becomes inaccurate, and therefore the CSI cannot be promptly and successfully used 

>. 

_ for relay selection protocols (RSPs). Correspondingly, the system performance drastically deteriorates. However, the 

• shadowing side information (SSI) varies more slowly such that its coherent time is substantially greater than both 

that of the CSI and backhaul /transmission delays in feedback channels. Within this concept to circumvent these 
. undesirable limitations, we suggest in this paper that some performance gain can be still obtained utilizing only 

the SSI for RSPs. In contrast to the other RSPs available in the literature, we propose a partial-RSP in dual-hop 
amplify-and-forward relaying system utilizing only the SSI of the relay terminals instead of their CSI. Furthermore, 

X 

^ ■ we offer a unified and generic performance expression of the proposed partial-RSP, spontaneously combining the 

analyzes of the average bit error probability, ergodic capacity, and moment generating function as a single analysis. 
Finally, as an illustration of the mathematical formalism, some numerical and simulation results are carried out for 
the generalized composite fading environments whose fading envelope follows extended generalized-K distribution, 
and these numerical and simulation results are shown to be in perfect agreement. 
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I. Introduction 

RELAY technology has recently gained great interest amongst researchers, engineers and developers 
in both academic and industrial fields as a transmission technique promising in millimeter wave (60 
GHz or above) radio frequencies, improving the quality of service, providing high data rate, and extending 
the coverage area without additional transmit power in deeply shadowed (composite) fading environments 
|[Il-|l3l. More specifically, referring to the current research and development projections we are currently 
dealing with, it is reasonably estimated that in future wireless technologies (5G and the beyond), there will 
be explicitly many wireless high-speed computers / devices / phones using external central process unit 
(CPU), external memory and external storage trough the medium of wireless communications links while 
their sizes are getting smaller and smaller. These applications naturally require high wireless data transfer 
with sophisticated usage of millimeter wave radio frequencies. Due to severe propagation characteristics 
in millimeter wave radio frequencies, the coverage area is generally within a 10 m or much smaller radius. 
Accordingly, wireless communications systems require to be assisted by relay technology such that the 
transmission between the source terminal (ST) and the destination terminal (DT) is successfully achieved 
by one or more relay terminals (RTs). There exist two famous relaying techniques one of which is the 
amplify-and-forward (AF) technique [3] in which each RT simply amplifies the received signal and then 
forwards the amplified signal to the DT. The other one is the decode- and-forward (DF) technique HI 
in which each RT decodes the received signal then forwards the decoded signal to the DT. While most 
frequently using these two relaying techniques, the purpose of relay technology is not only to increase 
the quality in signal transmission by obtaining spatial diversity at the DT via mutually independent copy 
of the transmitted information [[Il-[l3l, |l5l but also to increase the reliability of the signal transmission 
by reducing the fluctuations in the signalling throughput |l6l. For that reason, more than one RT can 
be employed in the signal transmission iTTl- lfTOl , but this solution considerably increases the network 
complexity with respect to the channel estimation, signal detection and transmission synchronization as 
well as information combining, resource allocation, scheduling and power control for signal transmission 
[fm . In addition, as the number of RTs increases, the signal processing load at the central entity (CE) 
and DT drastically increases [11]. In order to maintain the corresponding network complexity as simple 
as possible while increasing the overall quality of service and keeping the acceptable reliability, the 
best known solution is to select only one RT which is in the best fading conditions to help the signal 
transmission from ST to the DT d, [fTT1l-[fT3]|. 
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Designing a wireless communications system based on a relay technology considerably depends on 
which relay selection protocol (RSP) is utilized. Specifically, the various types of RSPs are well studied 
and discussed in the literature for different fading environments iTTll- inn . [[T3l - [|2n . These protocols can be 
grouped as best-relay selection Q-BUl, lfT3l . lfT4ll . IfTSl . best-worse-channel selection Q, best-harmonic- 
mean selection |I71, [[T4ll . nearest- neighbor selection ifTTI . ifTSl . reactive-relay selection [[T9ll . and partial 
relay selection lfT6ll . In more details, the best-relay selection suggests that the RT that can achieve the 
highest signal-to-noise ratio (SNR) at the DT is chosen from all the available RTs. For instance, lfT4l 
studies the error probability performance of this protocol, and an upper bound is presented in a closed- 
form for Rayleigh fading environments. The other RSP is the best- worse-channel selection [[TOl . IfTTI in 
which the RT with the largest of the worst channels from either ST — t- RT or RT — DT transmission links 
is selected among all the available RTs. The outage performance analysis for this protocol is successfully 
carried in both IfTTI and IfTOl for Rayleigh fading channels. It is further worth mentioning that both best- 
relay selection and best-worse-channel selection protocols achieve full diversity order. The other RSP 
achieving full diversity order is the best-harmonic-mean selection [fTl, lfT4fl in which the RT that has the 
largest harmonic-mean of the SNRs of the ST — t- RT and RT — DT fading channels is selected. It is 
important to note that in order to implement these RSPs, there needs to be a CE collecting all channel state 
information (CSI) of all available transmission channels from ST — > RTs and RTs — DT. Specifically, 
as the number of the RTs increases, the CE is drastically overloaded with all CSI information and then 
these RSPs, which maintain full-diversity order, become unfortunately impractical. 

In order to simplify the selection complexity and reduce the overload at the CE, the nearest-neighbor 
selection is proposed in IfTTI . IfTSl in which the RT nearest to the DT is selected, where "the nearest RT" 
is not necessarily the spatially nearest RT to the ST or DT but the RT with the strongest transmission 
channel (i.e., with the transmission channel having less amount-of-fading and less shadowing /path-loss) 
to the RT or DT. In addition to the nearest-neighbor selection, the other complexity-reduced RSPs such 
as reactive-RSP and partial-RSP were proposed in lfT6l . ifTOl - lflTI . respectively. Note that both reactive- 
RSP and partial-RSP require the CSI information of only one of channels ST — RT and RT — )• DT, 
respectively. Hence, the complexity of the network is substantially reduced by these two protocols. More 
specifically, the reactive selection is an quite efficient solution especially when RTs are close to the ST, 
and in which the RT is selected achieving the highest instantaneous SNR at the RT — t- DT channel |fT9l . 
On the other hand, the partial-RSP fflSl, |f20l, fflll in which the RT with the highest SNR at the ST RT 
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channel is selected turns into an efficient solution when RTs are close to the DT. 

As mentioned above, there exists a need in all RSPs to send all required CSI information to a CE 
through feedback channels in order to successfully achieve relay selection. However in practical cases 
in millimeter wave radio frequencies, there are some certain limitations for sending all required CSI 
information. For example, 

i) The number of available feedback channels in a communication system is typically limited and 
generally smaller than the number of RTs, which inherently turns into a feedback channel assignment 
problem specially when the coherent time of the CSI information is substantially small, 

ii) As the coherent time of the CSI information gets smaller, the energy consumption in the feedback 
channels significantly grows due to the throughput quality requirements. 

iii) The feedback channels, regardless of being wireless or wireline, have a backhaul /transmission delay 
which is practically so much higher than the coherent time of the CSI in the millimeter wave radio 
frequencies. As such in practical applications, the CSI will be inherently altered in the fading 
channel before successfully reaches to the CE through feedback channels. Specifically meaning 
that under these specified conditions, the CSI in the fading channel and the CSI reached at the 
CE will be assuredly uncorrelated. Thus, the CSI received by the CE through feedback channels 
explicitly becomes inaccurate for the RSPs, and consequently the system performance drastically and 
distinctively deteriorates. 

It is particularly useful to emphasize that these limitations even get worse as the mobility of the ST, RTs 
and DT increases. For example, increasing the carrier frequency in wireless systems proportionally and 
apparently magnifies the Doppler effects ||22|| . Following that the Doppler frequency in 60 GHz is typically 
at least 10 times greater than that in 6 GHz frequencies under the same mobility conditions [|22ll - [|24ll . 
which clearly states that there exists a proportional relation between Doppler's effect and the time-varying 
nature of the fading environments. Thus, the rate of the CSI information in 60 GHz frequencies is certainly 
at least 10 times greater than that in 6 GHz frequencies ||25]| . In other words, the CSI varies very fast so 
its coherent time is much smaller than that of a CSI in low-frequencies (about 6 GHz or below) under the 
same mobility conditions. It is therefore not difficult to envision such a practical scenario in millimeter 
wave radio frequencies that all RTs will continuously try to use at least one feedback channels but some 
RTs could not be assigned to feedback channels because of the limited number of the feedback channels. 
In such a case, the RTs have to reduce their CSI feedback to a low rate in order to share the feedback 
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channels amongst themselves. Because of reducing the feedback of the CSI, the CE will not correctly 
achieve the RSP, and the overall performance will take a turn for the worse. In particular under the 
circumstances when the CSIs of the RTs are non-identically distributed, the performance of the employed 
RSP is bounded by that of the average power (AP)-based RSP (i.e., pathloss-based RSP) in which the 
RT having the highest AP/the lowest pathloss is chosen from all RTS [fTTIl . Otherwise when CSIs are 
identically distributed, it is bounded by that of the round-robin (RR)-based RSP in which each RT is 
periodically chosen for signal transmission. 

In addition to the feedback channel assignment problem, the feedback channels, regardless of being 
wireless or wireline, have a backhaul /transmission delay which is practically so much higher than the 
coherent time of the CSI in the millimeter wave radio frequencies. As such in practical scenarios, the 
CSI will be inherently changed in the fading environment before successfully reaches to the central entity 
(CE) through feedback channels. Under these circumstances, even if there are enough number of feedback 
channels, the CSI in the fading channel and the CSI reached at the CE will not be same, and even within 
this practical scenario, they will certainly not be correlated. Consequently, the CSI cannot be promptly 
and successfully used for RSPs and the system performance explicitly and drastically deteriorates. Within 
this context, it is worthy noting that the shadowing side information (SSI) of a fading channel varies 
very slowly when compared with the CSI such that its coherent time is substantially greater than that 
of the CSI and backhaul /transmission delays in feedback channels. Specifically meaning that even if 
the CSI is surely difficult to be utilized for a RSP in millimeter wave radio frequencies, the SSI of the 
fading channel could be at least easily used to obtain some performance gain. As such, the relay selection 
based on the SSI information appropriately turns into an interesting solution for wireless communications 
systems operating in wireless millimeter wave frequencies. In particular, in contrast to the available RSPs 
ir/l- lfTTTl . [[T3l - [|2ni . we propose in this paper a new RSP based on the SSI instead of CSI, and we 
analyze its performance measures by offering a unified and generic performance expression that shows 
how to treat different performance analyzes as a single performance analysis. Explicitly, the followed 
analysis spontaneously combines the average bit error probability (ABEP), ergodic capacity (EC), and 
moment generating function (MGF) for a variety of fading environments. In this context, we reintroduce 
the extended generalized-K (EGK) distribution Il26ll , [|27l as a unified model for shadowed fading channels 
and then give some computationally efficient analytical results for performance measures. 

The remainder of the paper is organized as follows. In Section II, the system model of the partial- 
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RSP in dual-hop AF relaying system that only employs the SSI is overall described in details for 
generalized composite fading environments,Then, the unified and generic performance analysis is presented 
in Section III such that some main results regarding the special cases of the proposed unified and generic 
performance expression are briefly discussed. For possible applications in composite fading environments, 
the EGK distribution is reintroduced in Section IV as a generalized composite distribution that provides a 
unified theory to model the envelope statistics of wireless communication channels commonly used in the 
literature. In the following, its cumulative distribution function (CDF) and reciprocal moment generating 
function (MGF) are obtained in a closed form. Then for EGK fading environments, the correctness and 
accuracy of the unified performance analysis have been numerically checked. Finally, the conclusions are 
drawn in the last section. 

II. Partial Relay Selection Based on Shadowing Side Information 

We consider an AF dual-hop transmission in millimeter wave radio frequencies, between one ST S and 
one DT D through the L > 1 number of RTs {Ri, R2, . . . , Rl} as seen in Fig. [T] Assumed that there 
does not exist a direct link between the ST S and DT D due to the unsatisfactory channel quality, then 
the transmission can be performed only via RTs. In addition, appropriately assumed that the L RTs are 
relatively close to the destination. However as mentioned before, due to the impediment of the limited 
number of feedback channels, the backhaul / transmission delays subsisting in feedback channels, and the 
difficulties in synchronization and power control among RTs, only one RT having the best conditions is 
selected to enable the transmission between the ST and DT. More precisely as seen in Fig. [TJ for all 
hops A; G {1, 2} and RTs £ E {1,2,..., L}, {7^/} denote the composite fading coefficients defined as the 
product of shadowing {Sk/} and multipath fast fading {Sk/} such that 'jk/ = <Sk,iQk,e- Without loss of 
generality, the power of multipath fast fading, i.e., E[^fc £] is assumed to be unity, where E[-] denotes the 
expectation operator. Hence, the short-term and long-term average powers of the composite fading 7^. £ 
are IE[7fe^£] = Sk/ and ]E[7fc £] = E[5fc,£] = ^lk,£, respectively. Accordingly within this practical scenario, 
a RT could be in general picked from a set of RTs for the best transmission, archetypically by means of 
utilizing either CSI composed of shadowing and multipath fast fading components. 

However previously mentioned some certain limitations are, we propose a partial-RSP based upon the 
SSI comprising only shadowing component for communications systems operating in millimeter wave 
radio frequencies. Explicitly within this proposed partial-RSP, the ST S broadcasts predetermined training 
symbols to all RTs i?2, • • • , Rl}- Then, for all £ e {1,2,..., L}, the M RT estimates the SSI Si^e 
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with the aid of training symbols recovered from the transmission channel from the ST S. As a sequel 
using the feedback channels, each RT sends its SSI information to the CE which selects the RT having 
the highest SSI as 

M = arg max Si i (1) 

ee{i,2,...,L} ' 

where M E {1,2,. is the index of the selected RT. It is important to note that the index of the 

selected RT, i.e., M is a discrete random variable (RV) because of random nature of the SSI variation, 
and its probability mass function (PMF) can be simply written as 

L 

Pm (ni) = ^ /i£ 5 (m - , (2) 

where S (■) is the Dirac's delta function [|28l Eq. (1.8.1)], and where the coefficient /i£ denotes the probability 
of that the £th RT is selected, and it is namely defined as = Pr [M = £], i.e., 

^^e= PSi.i (s) Yl Psi.fc (s) ds, (3) 
Jo 



fc=i 



such that fJ-e = 1, where for s E M"*", both ps^ ^ (s) and Ps^ ^ (s) respectively denote the probability 
density function (PDF) and CDF of the SSI 5i ^ that the first-hop transmission is subjected to. Under the 
condition of that the Mth RT is selected, the ST S transmits the signal with a unit average power and 
the Mth RT receives the signal and amplifies it with the gain 

1 1 

Gm = = I (4) 

where Si^m — maxfcgji 2,...,l} Si^k denotes the shadowing of the first-hop fading channel conditioned on 
that the Mth RT is selected. In this sense, the instantaneous SNR of the first-hop transmission becomes 
7i,M = Si^M GiM- Then, remembering that the index of the selected RT M is a discrete RV with the PMF 
given in the overall instantaneous SNR 7^^^ at the DT D can be written in terms of a hyper-distribution, 
that is L ^ 

OM,e, (5) 

11/ + 12/ 

where 71 ^ = Si^cQi^e denotes the instantaneous SNR of the first-hop transmission conditioned on that 
the £th RT is selected, and where is the Kronecker's delta [28, Eq. (5. 10. 1.3)] such that 5„ ^ = 1 if 
n = k and 5n,k = otherwise. As mentioned before, the RTs are quite close to the DT D, creating a relay 
cluster. In this case, it is appropriate to assume that the transmission channels from RTs Ri, R2, ■ ■ ■ , Rl 
to the DT D are subject to negligible shadowing effect, which means that for all i E {1,2,..., L}, the 



^ ' Q ^ (>M,£- (6) 
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variance of the shadowing becomes Var [^2,^] ~ 0, where Var [•] is the variance operator, and its PDF 
P52f (■5) turns into a Dirac's delta distribution p^^^ (s) = 5 {s — i^2,£)- With this inference, the overall 
instantaneous SNR 7^^^ can be easily given as 

L 

^ 7l,£ + ^2,£G2,£ 

In the next sections, the performance of the SSI-based partial-RSP is analyzed, and some key results 
are presented with the discussion for generalized composite fading channels. 

III. Exact Performance Analysis over Generalized Fading Channels 

In this section, we study the exact performance analysis of the partial-RSP, depicted in Fig. [H using 
the SSI of all fading channels from the ST S to all RTs Ri, R2, . . . , Rl- 

A. Unified Performance Analysis 

The ABEP and EC are two important performance measures to quantify the achievable performance gain 
with the aid of the SSI-based partial-RSP. Over the past four decades, these ABEP and EC performance 
measures have been considered as two different problems. As such for these two performance measures, 
many separate and different solutions have been proposed in the literature for a variety of modulation 
schemes, diversity combining techniques in fading environments. The authors have recently offered in 
[|29l a generic unified performance (UP) expression that mathematically combines the ABEP and EC 
performance measures in a tractable way for a certain value of instantaneous SNR 7g„^, that is 



n I i-lTc 



'Pup hend) =C--<C- -fTTT^G 



l,n 



6,0 



(7) 



2 r v{h) 

where F (■) is the Gamma function [|30l Eq. (6.5.3)] and G^^"[-] denotes the Meijer's G function |I3T1 
Eq. (8.3.22)], and where the coefficient set Ai"\ n G N is defined as 



n times 



Ai") ^SX^'^To^. (8) 
It is informative to note that substituting n = 1 and c = 1 into Q results in the well-known conditional 
bit error probability (BEP) Pbep (lend) proposed by Wojnar in (321 Eq. (13)], that is 

Vbep{i.^,) = ^^, a,bell,l], (9) 



2T{b) ' ' t '2 

where F (-, ■) denotes the lower incomplete Gamma function Il33l Eqs. (8.350/2)], and where a depends on 
the type of modulation (^ for orthogonal frequency shift keying (FSK), 1 for antipodal phase shift keying 
(PSK)), b depends on the type of detection (| for coherent, 1 for non-coherent). Additionally noting that 
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for a fading channel with the bandwidth W Hz, substituting a = 1, b = 1, c = W/ log(2) and n = 2 into 
(|7]) where log (■) is the natural logarithm, (|7]) simplifies into the conditional channel capacity Pc (lend) 
bits/s/Hz, that is 

'Pc (lend) = W l0g2 (1 + 7e„J , bits/s/Hz, (10) 

where logg (■) is the binary logarithm. 

It is additionally useful to mention that the instantaneous BEP Pbep (lend) given by Q simplifies 
to Vbep (lend) — ^ i~ '^1 end) noH-cohcrent binary modulation schemes. With this inference, the 
unified expression Vup (lend) given by (|7]) can also be employed to obtain the MGF of the instantaneous 
SNR 7g„^. More specifically, setting a = p with p e C A ?R.{p} e R^, b = 1, c = 2 and n = 1 in ^ 
simplifies into the instantaneous MGF, i.e., 

PmGF hend) = ( -p7end) . ( 1 1 ) 

which enables the utilization of (|7]) so as to obtain the statistical characterizations such as PDF, CDF and 
moments within the same performance analysis framework. 

Eventually, using the unified performance expression Vup (lend) given in (|7]), the average unified 
performance (AUP) Vaup = ^[Pup (lend)] of partial-RSP based on the SSI of the first-hop fading 
channels available between the ST S and the RTs Ri, R2, ■ ■ ■ , Rl is given in the following theorem. 

Theorem 1 (Unified Performance). The AUP Vaup of the SSI-based partial-RSP is given by 

Kb,c iu)J2fiAMi/^,^, (u) -Q^^ihnx (^) + -^1/72,. («) ^^i/7i,£ («) ^^^^ 
where the auxiliary function Z^i, (u) is defined as 



Z 



n 



n J (-1) c 1 

[U] = C <C - ^ TT—G 

2] T{b) 



n'i 



au 



6,0,0 



(13) 



where a, b, c and n are the performance parameters explained before. Furthermore in (|T2]) . for all 
i E {1,2,..., L}, the coefficient fi£ is given in ^ denoting the probability of that the ith RT is selected. 
Further, both Mi/^^ (p) = ]E[exp (— p/71 £)] and M-ih^i (p) ~ ^^[^xp (— J9/72/)] denote the reciprocal 
MGF s for the first and second hops' instantaneous SNRs 71,^ and 72,^. 

Proof: Explicitly for a certain non-negative distribution of the instantaneous SNR 7e„^, the AUP 

Paup = ^[Pup (lend)] is Completely represented as 



n 



'^--- = '-2V+ Tib) .0 



p l,n 



aS"\i 

6,0 



-^My (p) 

Qjj 'end ^'^ ^ 



dp 



(14) 
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using the MGF-approach given in ||29l Eq. (30)]. Furthermore in (fT4l) . Ai^ ^{p), p E denotes the 
MGF of the overall instantaneous SNR 7g„^, and is in general impossible to be in closed form because of 
the nature of composite fading channels. However in order to reduce (fT4l) into a computationally efficient 
single-integral, Ai^ ^(p) is required in a closed-form. Averaging with the PMF given by Q, the overall 
MGF M^^^aip) = E[e~-?^^endj can be written as 



M 



1=1 



exp —p 



7i,e'l2,£ 
7i,f + 12/ 



(15a) 
(15b) 



where for all £ G {1,2, .. . ,L}, the coefficient fii which denotes the probability of that the ith RT is 
selected is given in ([3]). Note that the instantaneous SNR at the DT becomes the normalized harmonic 
mean of the first and second hops' instantaneous SNRs 71^^ and 72,^ when the £th RT is selected. More 
precisely, after performing some algebraic manipulations using ll34l Theorem 3], Ai-y^^^ip) obtained as 

Jo (2v^) 5^/^^ I -^1/71,^ («) ^-^1/72,^ (w)+-^i/72,. («) ^-^V7i,f iu)\du, (16) 

where Jq (') denotes the zeroth-order Bessel function of the first kind [|33l Eq. (8.411/1)]. Eventually, 
substituting (fT6l) into (fT4l) and changing the order of the resultant integrals, and then performing algebraic 



manipulations using {8/ dp) Jq = —\fuJpJx {'^Vp^)' we have 



n 



Vaup = c - - <( c + 



(-iTc 

m 



/i^x {u) X 



i=i 



-^1/71,^ (^) ^-^1/72,. ('") + -^1/72,. («) ^-^1/71.^ (^) 1^^^ ^ (17) 



where the function X(m) can be obtained in closed-form by means of utilizing 1131] Eq. (2.24.4/1)] and 
[|3T] Eqs. (8.2.2/14) and (8.2.2/15)] as 



^0 VP 



a 


Ai"\l" 






aS")" 




aw 


P 


6,0 




6,0,0 



(18) 



Then, substituting (118]) into (1191) results in 



Vaup = c - - < c 



£=1 



an 



6,0,0 



X 



|-^l/7i.. («) ^-^1/72,^ («) + -^1/72,. («) ^-^1/71,. (19) 
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where the terms outside of the integration can be easily put inside the integration by using the well-known 
fact that roo r Q d ^ 

1^1/71,. (^) ^-^1/72.. (^) + -^1/72,. (^) ^^1/71.. («) = (20) 

Finally, applying (l20l) on (fT9l ) readily results in (fT2l) . which proves Theorem [TJ ■ 
It is both useful and informative to notice that the exact AUP Vaup presented in Theorem [T] is straight- 
forward to utilize depending on the exact closed-form expressions for the reciprocal MGFs Aii/^^i, (p) and 
A^2/72£ (p)- However, both 71 ^ and ^2/ are two composite distributions. Specifically, the instantaneous 
SNR 7i_f is given by 7i^£ = Si^iQi^e, then on the condition that the SSI knowledge, i.e., Si^e is available. 



[s)ds is written as 



the reciprocal MGF A^i/^,, (p) = E e~P/"''''\Si,i = s 

roo 

^i/ii.i (P) = / ^i/Gi.i ip/s)p§^^{s) ds, (21) 
^0 

where Aii/g^ ^ (p) is the reciprocal MGF of the multipath fast fading that the £th RT is exposed to. In turn 
as explained before. Si/ = maxke{i,2,...,L}<Si^k denotes the shadowing power of the first hop conditioned 
on that the £th RT is selected. Its PDF pg^ ^ (s) is given by 

L L 

P§.A') = Y.P^.A^)Up^.A')- (22) 

1=1 fe=i 

k^l 

where for £ G {1, 2, . . . , L} ^ (s) and P^^ ^ (s) are defined before. Note that referring |[27l , Mi/g^ ^ (p) 
results in closed-form expressions for the fading distributions commonly used in the literature. As a result 
of that, the derivative of A^i/g^ ^ (p) is also available in closed-form in the literature. Thus, the complexity 
on the computation of the reciprocal MGF Aii/j^,, (p) is simplified by all means. As a step forward, using 
Gauss-Chebyshev quadrature (GCQ) formula QUI Eq.(25.4.39)], the reciprocal MGF A^i/^^^ (p) can be 

written as ^ 

^1/11.1 (P) = Vn,eMi/g, , {p/ Sn) , (23) 

n=l 

where the coefficient r]n,e is given by 

L L 

'r]n,£ = Wn 

/ = ! k=l 

Furthermore in (l23l) and (l24l) . the coefficients s„ and Wn are defined in [[35l Eqs. (22) and (23)], respectively, 
and where the truncation index N could be sufficiently chosen for an enough accurate result. In addition, 
conditioned on that the £th RT is selected, the instantaneous SNR of the second hop, i.e., 72/ is given by 
72,£ = S2/Q2,t- Then, the reciprocal MGF Mi/^^^ (p) = E[exp (— ^5/72,^)], ^{p} G IR+ is obtained as 
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-^1/72.. (P) = / 02,1 (p/s) PS2,, (s) ds (25) 

^0 

where A^2/g2 e (p) reciprocal MGF of the muhipath fast fading that the ith RT is exposed to. In the 

special case of that the RTs R2, R2, ■ ■ ■ , Rl are much more closer to the DT, which create a cluster of 
RTs whose channels towards to the DT D have very negligible power fluctuation as compared with their 
average powers. In other words, Var [^2,^] ~ for all RTs i E {1,2,..., L}. Therefore, the distribution of 
the shadowing S2/ turns into Dirac's distribution and its corresponding PDF ps^ ^ (s) definitely becomes 
Ps2e i^) = ^(■^ ~ ^2/)- Immediately substituting this PDF into (|25] ). -Mi/^^e (p) results in 

-^1/72,. (P) = -^i/G2.. {p/^2,i) ■ (26) 

Finally, note that (l23l) . (|24|) . (|25l) and (|26|) can be efficiently computed for a variety of commonly 
used fading distributions such as Rayleigh, lognormal, WeibuU, Nakagami-m, generalized Nakagami-m, 
generalized-K and EGK distributions. Consequently, substituting them into (fT3] ). the exact AUP Vaup is 
easily and efficiently computed. 

B. Special Cases of the Auxiliary Function 2^"^^ (u) 

Let us consider in this section the special cases of the auxiliary function Z^'^ ^ (n) in order to check 
and show analytical simplicity, correctness and accuracy: 

Special Case 1 (Auxiliary Functions for the Bit Error Probabilities of the Coherent Modulation Schemes). 
Setting n = 1, c = 1 and b = 1/2 in (fT2l) in order to find the auxiliary functions for the bit error 
probabilities of coherent BFSK and BPSK modulation schemes, the auxiliary function Z^'^ (u) simplifies 
to 

^a,i/2,i iu) = l- -Si (2v^) , (27) 

by means of using [|3T1 Eq. (8.4.12/1)], where Si(-) is the sine integral function defined as Si (z) = 
u~^sm{u) du [|30l Eq. (5.2.1)], and where the parameter a = 1/2 for orthogonal coherent FSK (BFSK) 
Il36l Eq. (8.43)], and a = 1 for antipodal coherent PSK (BPSK) [|36l Eq. (8.19)]. Inserting ^ into ([H]). 
the ABEP of the partial relay selection is obtained for BFSK and BPSK modulation schemes. □ 
Special Case 2 (Auxiliary Functions for the Bit Error Probabilities of the Non-Coherent Modulation 
Schemes). In the special case of non-coherent modulation schemes, substituting n = 1, c = \ and 6=1 
(fT3]) simplifies the auxiliary function by means of using ||3T1 Eq. (8.4.19/1)], that is 

^ii,iH = ^^(2v^), (28) 
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where the parameter a = 1/2 for orthogonal non-coherent FSK (NCFSK) Il36l Eq. (8.69)] and a = 1 for 
antipodal differentially coherent BPSK (BDPSK) |l36l Eq. (8.85)]. Eventually by substituting ([281) into 
(fT2l) . the ABEP of the partial relay selection is easily obtained for non-coherent NCFSK and BDPSK 
modulation schemes. □ 
Special Case 3 (Auxiliary Functions for Bit Error Probabilities of the Wojnar's Representation). For 
conditional BEP Pbep (lend) proposed by Wojnar in [[32l Eq. (13)], setting n = 1 and c = 1 in (fT2l) . the 
auxiliary function expressly simplifies to 

1 1 



7^ 



U 



2m 



^1,3 



au 



1 

6,0,0 



(29) 



After performing some algebraic manipulations either by means of using OTl Property 2.11] or by 
recognizing that the CDF of a Fox's H distribution can be written in two different forms given in ||38l 
Eq. (4.17)] and JMl Eq. (4.19)], ^ simplifies to 

1 



7^ 



U 



2T{h) 



^-^1,3 



au 



1 

6,0,0 



(30) 



which is in perfect agrement with the well-known result given in [1391 Eq. (14)]. In more details using Oil 
Eq. (8.4.12/1) and (8.4.19/1)], substituting 6 = 1/2 and 6 = 1 into ^ result in ^ and (EHl), respectively. 
Finally, inserting (|30l ) into (IT2l) . the ABEP of the partial relay selection can be readily obtained for all 
coherent and non-coherent binary modulation schemes. 

Special Case 4 (Auxiliary Function for Ergodic Capacity). Substituting a = 1, 6 = 1, c = W/ log(2) and 
n = 2 into ([T2l) . the auxiliary function for the EC analysis is easily obtained as 



1,1, 



w 

log(2) 



[U] 



log(2) 



2,3 





1,1 " 


u 






1,0, 0_ 



(31) 



Note that utilizing [HOl Eq. (07.34.03.0475.01)] and [|40l Eq. (06.06.03.0003.01)] together, (HB further 
simplifies into [1411 Eq. (20)], that is 

W 



1 '■ 1 /ON 



(log(M) — Ei 



-u 



C) 



(32) 



where Ei (■) is the exponential integral function ll33l Eq. (8.211/1)] and C is the Euler-Mascheroni constant 
(also called Euler's constant) |l33l Eq. (8. 367/1)]. Eventually, inserting (|32l) into (fT2l) . the EC of the partial 
relay selection can be readily obtained. □ 
Special Case 5 (Auxiliary Function for the Moments-Generating Function). Substituting a = p, h = \, 
c = 2 and n = 1 into (fT2l) . the auxiliary function for the MGF analysis is easily obtained as 
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^1,1 
"^1.3 



pu 



1 

1,0,0 



(33) 



Utilizing gOl Eq. (07.34.03.0285.01)], ^ Eq. (7.14.2/84)] and [gOl Eq. (9.6.3)], (133]) further simplifies 
into the transform function of [|42l Eq. (13)], that is 

^ii,2 (^) = ^0 (2v^) • (34) 
Finally substituting (l34l) into (fT2)) yields (fT6l) as expected. □ 
Consequently, the mathematical formalism in this subsection shows how to treat the ABEP analysis of 
coherent and non-coherent binary modulation schemes, EC analysis and MGF-based statistical analysis 
of the overall instantaneous SNR simultaneously in parallel. As such, these different analyzes could be 
evidently viewed as a single problem ||29l for a variety of fading distributions. 



C. The Moments of the Overall Instantaneous SNR 

It is worthy to mention that the moments of the overall instantaneous SNR 7^,^^ is another important 
tool to characterize the probability distribution such that it can be used to obtain the higher-order amount- 
of-fading (AF) defined as [|43l Eq. (3)] 



endi 



(35) 



end\ 



where the A;th moment of the overall instantaneous SNR 7^^^ can be obtained with the aid of ]E[7g, 



end} 



{d/dpYM..^ ^{p)\p^Q^ where M.^ ^{p) is given in (fT6l) . Within this concept, (|34|) is very useful. Ac- 
cordingly, the kth derivative of (|34|) with respect to p at p = is derived by means of performing some 
algebraic manipulations using [[331 Eq. (8.411/1) and (8.411/2)], that is 



[u 



p^O 



u 

Id 



(36) 



Utilizing (l36l) . the kth moment of the overall instantaneous SNR 7g„^ is easily obtained as 



niLd] 



d 



Then with aid of (137]) . the higher-order AF can be readily obtained for a variety of fading distributions. 
In addition, (l37l) can be utilized in either asymptotic performance analysis [[431 , [[44l or approximation of 
the outage probability using Laguerre moments [[34l . 

Despite the fact that the novel and generic technique represented in this section is easy to use for 
a variety of generalized fading distributions, let us consider in the next section some numerical and 
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simulation examples achieved for the EGK fading environments in order to check its analytical simplicity 
and accuracy. 

IV. Applications in Generalized Composite Fading Environments 

We consider the partial-RSP based on the SSI of the fading channels between the ST S and the 
RTs Ri, R2, ■ ■ ■ , Rl, as depicted in Fig. [T] for composite fading environments. However, even if the 
results obtained in the previous sections are in general applicable to many generalized composite fading 
environments, our results in this section have been intentionally carried for the EGK fading distribution 
in order to check numerically the correctness and accuracy of the presented mathematical formalism. In 
the following subsection, the EGK distribution model is explained and its some important properties are 
accentuated. 

A. Extended Generalized-K Fading Environments 

The EGK distribution is a generalized distribution that provides a unified model to statistically character- 
ize the envelope statistics of wireless communication channels as explained in [|26l . [|27l . More precisely, 
as listed in both ll26l Table I] and [f27l Table I], several distributions such as Rayleigh, lognormal, WeibuU, 
Nakagami-m, generalized Nakagami-m, generalized-K are special or limiting cases of the EGK distribution. 
Referring to the partial-RSP depicted in Fig. [U for all hops k e {1,2} and RTs £ e {1, 2, . . . , L}, the 
instantaneous SNR 7^. ^ is assumed to follow the EGK distribution whose PDF can be given in [|26l Eq.(3)], 
m Eq.(5)] by 



where m^,^ (0.5 < m^,^ < 00) and .^^.^ (0 < £ < cxd) represent the fading figure (diversity severity 
/ order) and the fading shaping factor, respectively, while ra^,^ (0.5 < < 00) and C^,^ (0 < Ck,e < 
00) represent the shadowing severity and the shadowing shaping factor (inhomogeneity), respectively. 
In addition, the parameters 0;.^ and cpf^^^ are defined as 0;.^ = r(m;.^ + l/^j,^)/r(m^,^) and (ff^^^ = 
r('^fc£ + l/Cfc£)/r(n;.^), respectively. Furthermore, in (l38l) . T (-, -, -, ■) is the extended incomplete Gamma 
function defined as T (a, x, b, (3) = r""^ exp (— r — 6r~^) dr, where a, /3,b E C and x E M"*" Il45l Eq. 
(6.2)]. 

In addition, its MGF, i.e., M^,^^{p) = E [exp(-p7fc^)] is given in |[271 Eq. (15)]. Then, using ll34l 
Theorem 1], its reciprocal MGF Mi/j^^{p) = E [exp(— p/7^,^)] is obtained as 
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1 --3,0 



~n — P 



(0, 1), (m.^, l/^fc,£), {Uf^p, 1/Ck,£) 



(39) 



where H;;;^[-] denotes the Fox's H functiorUBI defined in ^ Eq.(8.3.1/1)], JSTl Eq.(l.l.l)], and where 

means that the coefficients are absent. 

As mentioned before, the instantaneous SNR 7^ ^ is defined as jk/ = <Sk/ Gk,e- When jk/ follows the 
EGK distribution, then its multipath fast fading part, i.e., Qk/ follows the well-known generalized gamma 
(GG) distribution with the PDF pg^, ^ (g) = E[5 (g — ^k,e) \ <Sk,e = 9 ^ where assumed without loss 
of generality that the multipath fast fading Qk,£ has unit power, i.e., E[^a,s^] = 1. Then, the PDF pg^^ (g), 
g E is explicitly given by 

P.,. (9) = |^(0.,.^)--^^--^e-(^--)^^^ (40) 

Note that as indicated in Il39ll , the reciprocal MGF ■M.i/g^^{p) = E,[exp {—p/Qk,e)] and its first-order 
derivative {d / dp)M.i/g^ (p) = —E,[exp{—p/Qi:,e) /Gk/] are required to obtain computationally efficient 
analytical results. With the aid of the result in ll39l Eq. (18)], the reciprocal MGF of the multipath fast 
fading Q^^^ is written as 

-^i/Sfc,, (P) = p/ ^ {m^p 0, ^k/P, l/^k,i) ■ (41) 
y^ki) 

Accordingly using Il45l Eq.(6.22)] and [[3T1 Eq. (8. 2. 2/30)] together, its first-order derivative can be written 
as (391 Eq. (19)] 

^^VSu, iP) = -f^f K - Ve., 0, <Pk,iP, 1/4,.) • (42) 

Note that the average power of 7^^, i.e., lE[7fc_^] = S^/ changes slowly, so without loss of generality 
Sk,e can be assumed a constant for a short period of time but has a distribution over a long period 
of time. Explicitly, its distribution explicitly follows the PDF given by p^^.^ (s) = E[5 (s — E[7fc^^])] = 

E[6 {s - Sk,i)], given by 



^("^kl) ^kl \^kl 



PS,, is) = X e V".. (43) 



where E[7fc/] = E[5fc/] = fi^, ^ for a long period of time. As mentioned in the previous subsection that 
the selection of the RT will be achieved based on the SSI information of the channels / hops. The CDF 

of Sk^e, i.e., Ps^^^is) = E[6 (s — E[7fc^£])] = E[6 {s — 3^/)], s G is therefore required and given by 

'Using 1311 Eq.(8. 3.2/22)], the Fox's H function can be represented in terms of tlie Meijer's G function 1311 Eq.(8.2.1/1)] which is a 
built-in function in the most popular mathematical software packages such as Mathematica® 

^The authors gives in the appendix of 1461 an efficient implementation of the Fox's H function in MATHEMATICA® 
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Ps,, {s) = -1^7 k,., {s , (44) 

where 7 (-, ■) is the upper incomplete Gamma function [[33l Eqs. (8.350/1)]. 

B. Unified Performance Analysis in Extended Generalized-K Fading Environments 

Regarding the numerical and simulation results, we considered the EGK fading environment for different 
number of RTs, i.e., L G {1,2,3,5} with the statistical channel parameters m^^, C^^, n^^^, and VLj^^^ 
for all hops k G {1,2} and all RTs £ G {1, 2, . . . , L} as shown in Table HI In order to obtain the AUP 
Vaup of the SSI-based partial relay selection in EGK fading environments, Theorem [T] introduces an 
MGF-based approach that requires only the reciprocal MGFs Mi/-,-^^ (p) and A^i/72f (p) in closed-form. 
As such, the reciprocal MGF Mi/^-^^^ (p) can be obtained by substituting (|42|) and (|44|) into (l23l) . that is 

N 

V,.,, iP) = E f^r(™,,„ 0. ±) . (45) 

where the coefficients {r]n/} are easily computed substituting both (l43l) and (|44)) into (|24)) . Accordingly 



using (1421) . its derivative can be easily given as 

N 



In addition, conditioned on that the ith RT has been selected, the reciprocal MGF of the second hop, 
Mi/y2e (p) '^^^ be obtained using (l39l) , that is 



1 --3,0 



02,^ V52,£ 
^ 



(0, 1), (m2,, 1/6,^), (^^2,^' VC2,£) 



(47) 



As mentioned before, the RTs Ri, R2, . . . , Rl are, without loss of generality, assumed to be very close 
to the DT D, accordingly creating such a RT cluster that the transmission channels from RTs to the DT 
D are not subjected to disruptive shadowing effects as explained in Section ??. Evidently in such a case, 
the shadowing conditions certainly disappear since the shadowing figures {rig^} explicitly approach to 
infinity when the RTs Ri, R2, . ■ ■ , Rl get closer to the DT D. Then refering (|4T)) . the reciprocal MGF of 
the second hop, M.i/^^^ (p) reduces to Mi/g^^{p/^2,e), that is 



Accordingly using (I42|) . its derivative can be easily given as 



7r-^V72,b = Vr;~^l"^2,£- ^.0,— — ,— . (49) 

dp T{m^f)n2/ V ^2.£ ^2/ ^2 
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Eventually, substituting the set of Equations ([45]), (|46l), (|48]) and (|49l) into ([HI), the AUP Vaup of the 
SSI-based partial-RSP in dual-hop AF system over EGK fading environment can be readily obtained as 

TV L 



— r(mi^- — ,0, ,— jrfm2_^,0,— — ,— 



+ 

—r( 772^^,0, , — ]T[m^^-—,0,— — ,— ] Uu. (50) 



Note that the coefficients /i^ in (|50l) . which is the probability that the ^th RT is selected, can be readily 
computed by substituting both (l43l) and (|44|) into ([3]). For the well-known generalized-K (GK) fading 
environment (i.e., for the special case of the EGK fading environment with the fading and shadowing 
shaping parameters = = 1 for all k e {1,2} and £ e {1,2,..., L}), the unified result given in 
(l50l) can be more simplified to 



r(mi^^)r(m2^^) ^ (l]2,^/m2j 



m2£ / m^^u\ I m^^u 



by means of using [|45l Eq. (6.41)], i.e., 

r (m, 0, 6, 1) = 2fe'"/2^^ (^2^6) , 3ft{6}GM+, (52) 

where K^^ (■) denotes the mth order modified Bessel function of the second kind [l30l Eq. (9.6.2)]. Note 
that, under the further assumptions that the shadowing distributions {Si^i} are identically distributed (i.e., 
n^ f^ = n, = C and Vli^i = Vl for all £ G {1, 2, . . . , L} such that Si^^ Si for all ^ G {1,2,..., L}), 
each RT will have the same probability to be selected, i.e., the probability of the £th RT is selected 
certainly becomes 



/if 



PsAs)[PsAs)Y~^ds = -, (53) 



for all ^ G {1,2,...,L}. In such circumstances, when the shadowing conditions get better (i.e., the 
disruptive effects of the shadowing conditions on the received power decreases) as the shadowing figure 
surely goes to infinity (n — oo). Then, the SSI-based partial-RSP turns into the simplest round-robin (RR)- 
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based partial relay selection tournaments in which each RT is periodically chosen for signal transmission. 
Although it has worse performance than the other RSPs, it provides quite an ABEP upper-bound (i.e., 
worst-case) performance that the SSI-based partial relay selection can achieve. In such a case, (l50l) 
simplifies as expected to 



^("^2,£>0, — — ) + 



P.„P = r zi, (u) ± i ^r(,n,, - -L. 0. -L)r( 

-r mi£,0,— — ,— r m2,^- — ,0,— — ,— ] Uu. (54) 
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On the other hand, when the shadowing distributions are not identically distributed (i.e., ^li,£ ^ 

^i,k, n^^i 7^ n^ ^, and Ci,£ 7^ Ci,fc for ^ 7^ ^^d when the shadowing figures goes to infinity (for all 
£ G {1,2,..., L}, ri^ ^ — 7- oo), the CE will definitely select the RT having the highest AP, that is 

I = arg max Qi g (55) 

£G{1,2,...,L} 

where i G {1,2, is the index of the selected RT. Note that the average power of each link is 
assumed to be constant and does not vary with data transmissions without loss of generality. Then, the 
index of the selected RT will not be a discrete random variable, and the SSI-based relay selection protocol 
will definitely turn into the well-known AP-based partial-RSP which is also known as pathloss-based relay 
selection |1 1| . In this case, (l50l) reduces into the unified performance of the AP-based partial-RSP, that is 




As a result from these analytical expressions, the performance of the SSI-based partial relay selection 
certainly deteriorates approaching to those of RR-based and AP-based partial-RSPs as the shadowing 
conditions in the fading channels from the ST to the RTs i?i,i?2, • • • , Rl get better. It is at this point 
worth mentioning that the shadowing conditions in all available channels from the ST to the RTs gets 
worse as the RTs get farther away from the ST. Then under severe shadowing conditions, the SSI-based 
partial relay selection turns into a feasible solution in the millimeter wave radio frequencies to maintain the 
performance of dual-hop AF relaying systems subjected to severe limitations such as fast the time-varying 
nature of the CSI and the backhaul / transmission delay in feedback channels. 
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C. Numerical Results in Extended Generalized-K Fading Environments 

As an illustration of the numerical accuracy of the presented mathematical formalism, we consider the 
EGK fading environments for different number of RTs, i.e., L E {1,2,3,5} with the statistical channel 
parameters m^^, n^^, ^^d for all hops k e {1,2} and all RTs £ G {1,2, ... ,L} as shown 
in Table H Accordingly, the ABEP performance of the SSI-based partial relay selection, obtained by 
simulation and numerical computations, is depicted in Fig. |2] and Fig. [3] for the BPSK and NCFSK 
modulations schemes, respectively, for different number of RTs. As seen in these figures, the numerical 
and simulation results are in perfect agreement. Further note that in our computer-based simulations, the 
number of samples for the instantaneous SNR 'jend is chosen 10^ in order to obtain quite accurate results 
for comparison reasons. 

With the fact that the performance of all RSPs significantly worsens to those of AP-based and RR-based 
RSPs because of some certain limitations mentioned previously, the performance of both the RR-based 
and AP-based partial-RSPs are also depicted in both Fig. [21 and Fig.[3l For example, the AP-based and RR- 
based partial-RSPs do not provide a significant performance gain while the average SNR and the number 
of RTs increase. On the other hand, the SSI-based partial-RSP has still good and respectable performance 
in shadowed fading environments. More specifically in Fig. [2] for a specific desired BPSK ABEP 2x10"^ 
achievable with L = 4 number of RTs, the SNRs of the AP-based and RR-based partial RSPs are 55 dB 
and 34 dB, respectively. On the other hand, the SSI-based partial-RSP reaches the same ABEP for only 
the average SNR 13dB. In addition in order to further quantify the performance gain of the SSI-based 
partial-RSP and demonstrate how it turns into a feasible solution in shadowed fading environments, we 
also compared in both Fig. [2] and Fig. [3] its ABEP performance with that of the CSTbased partial RSP 
even if the CSI-based partial-RSP is somehow difficult, as discussed in more details in the introduction, 
to realize in millimeter wave radio frequencies. Moreover for L = 4 number of RTs, we can see in Fig. [3] 
that the NCFSK performance of the SSI-based partial-RSP is approximately 5.5 dB worse than that of 
the CSTbased partial-RSP for the specific ABEP 10^'^. However, it provides some significant strengths 
to some certain limitations mentioned previously. Further with respect to the estimation, synchronization 
as well as information combining in millimeter wave frequencies, it maintains the network complexity 
as simple as possible while reducing the signal processing load at the CE and increasing the overall 
quality of service. Note that the other information theoretical framework for studying performance limits 
in the SSI-based partial-RSP is the channel capacity whose first-order statistics is typically known as the 
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EC denoting the maximum information rate for which arbitrarily small BEP can be achieved in fading 
environments. Having been obtained by SSI-based partial RSP, the improvement in the EC is illustrated 
in Fig. HI for different number of RTs, and compared to those of the RR-based and AP-based RSPs. As 
seen in Fig. |4l the EC of the SSI-based partial RSP improves as the number of RTs increases. 

In addition to this performance gain and complexity reduction, it is clearly expected that the performance 
gain of the SSI-based partial-RSP strictly depends on the random nature of the shadowing conditions. In 
particular, shadowing effects certainly disappear in the fading environments between the ST and RTs when 
the ST gets closer to the RTs. Under the same circumstances in EGK fading environment, the shadowing 
figures of the first-hop fading channels grow substantially (i.e., ^ — )■ oo for all i E {1,2,..., L}). For 
4 number of RTs under the reasonable assumption that their channels from the ST are exposed to the 
same shadowing distribution (i.e., n^^ = n for all i E {1,2,..., 4} and the other channel parameters 
are listed in Table U), Fig. [5] presents the ABEP of the SSI-based partial-RSP for the BDPSK modulation 
scheme with comparison to those of CSTbased, RR-based, and AP-based RSPs. Simulation and numerical 
curves show that the ABEP of the SSI-based partial-RSP gets gradually worsen to that of the AP-based 
partial-RSP as the shadowing severity in the first-hop fading channels decreases. Takin into account that 
performance behavior, we can quite easily deduce that the SSI-based partial-RSP certainly turns into a 
feasible solution in severely shadowed fading environments. 

V. Conclusion 

In conclusion, the RSPs available in the literature are typically using the CSI. On the other hand in the 
millimeter wave radio frequencies, the CSI varies very fast so its coherent time is much smaller than that of 
the CSI in low-frequencies under the same mobility conditions. Because of some certain limitations such 
as backhaul /transmission delays in feedback channels, the efficient usage of CSI for RT selection may 
not always be possible. However and fortunately, the SSI variates more slowly, so we showed in this paper 
that the SSI can also be employed in the partial-RSP to obtain some respectable performance in deeply 
shadowed fading environments. More specifically in order to achieve the RT selection in dual-hop AF 
relaying system, we proposed a partial-RSP utilizing only the SSI of the fading channels from the ST to 
the RTs instead of the CSI of these channels. Then for the computationally efficient numerical analysis of 
the proposed SSI-based partial-RSP, we offered a unified and generic UP expression which simultaneously 
combines the ABEP analysis of binary modulation schemes (BPSK, BDPSK, BFSK and NCFSK), EC 
analysis, and MGF-based statistical characterization in parallel. As such, this generic UP expression can 
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be efficiently computed for different performance analyzes and it is certainly an MGF-based approach 
applicable for a variety of fading environments. In addition, we reintroduced the EGK distribution for 
possible applications of the SSI-based RSP in composite fading environments. As an illustration of the 
mathematical formalism, some numerical and simulations are carried out and they are shown to be in 
perfect agreement. 
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Fig. 1. Partial relay selection with respect to shadowing side information. 
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Fig. 2. BPSK performance comparison in EGK fading environments. 
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Fig. 3. NCFSK performance comparison in EGK fading environments. 
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Fig. 4. Ergodic capacity comparison in EGK fading environment. 
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Fig. 5. For the average SNR 30 dB and 4 number of RTs, the BDPSK performance comparison with respect to shadowing figure of the 
first-hop fading channels in EGK fading environment (i.e., tij £ = n for all I £ {1,2,..., 4}). 



TABLE I 

Statistical Channel Parameters in Double-Hop Communications Operating in EGK Fading Environment 
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